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Much of the progress in our understanding of dynamo mechanisms has been made within the 
theoretical framework of magnetohydrodynamics (MHD). However, for sufficiently diffuse media, the 
Hall effect eventually becomes non-negligible. We present results from three dimensional simulations 
of the Hall-MHD equations subjected to random non-helical forcing. We study the role of the Hall 
effect in the dynamo efficiency for different values of the Hall parameter, using a pseudospectral 
code to achieve exponentially fast convergence. We also study energy transfer rates among spatial 
scales to determine the relative importance of the various nonlinear effects in the dynamo process 
and in the energy cascade. The Hall effect produces a reduction of the direct energy cascade at 
scales larger than the Hall scale, and therefore leads to smaller energy dissipation rates. Finally, we 
present results stemming from simulations at large magnetic Prandtl numbers, which is the relevant 
regime in hot and diffuse media such a the interstellar medium. 
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I. INTRODUCTION 

The generation of magnetic fields by dynamo activity 
plays an important role in a wide range of astrophysical 
objects, ranging from stars to clusters of galaxies. The 
gas in these objects is characterized by turbulent flows, as 
shown for instance by scintillation observations of the in- 
terstellar medium [l|, y] , or from pressure maps in galaxy 
clusters 0] . Mechanisms able to generate magnetic fields 
by dynamo action are often classified as large- and small- 
scale dynamos, depending on the correlation length of the 
induced magnetic field. In this context, large and small 
are referred to the energy containing scale of the turbu- 
lent hydrodynamic fiow. This classification is not rigid, 
as in many astrophysical objects both dynamos may be 
at work, but it gives a useful framework considering the 
limitations in the scale separation that can be achieved 
in numerical simulations. Also, the physical properties 
of the fiows that can give rise to one or the other are 
somewhat different. 

Helical fiows have proved efficient in generating large- 
scale dynamos, i.e., on scales larger than the energy- 
containing eddies of the flow [3t3|- It is now known 
that large-scale dynamo action can also be produced by 
anisotropic and inhomogeneous flows (e.g., flows with 
a large scale shear). On the other hand, non- helical 
flows can be instrumental in generating small-scale dy- 
namos Q, i.e., on sizes smaller than those of the energy- 
containing eddies pj|-ill|. In recent years, the study 
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of small-scale dynamos with magnetic Prandtl number 
Pm ~ v/rj (the ratio between the viscosity and the mag- 
netic diffusivity of the plasma) different from unity has 
received special attention [ij, dM i both for Pm ^ 1 [IJ] 
and for Pm <C 1 [ll, [3 ■ Motivations to study these 
regimes include recent experiments of dynamo action us- 
ing liquid sodium [l7| , as well as the fact that many astro- 
physical plasmas are characterized by magnetic Prandtl 
numbers different from unity. For instance, the magnetic 
Prandtl number is much smaller than one in the solar 
convective region, and it is typically much larger than 
one in the interplanetary medium and also in the inter- 
stellar medium (ISM). 

For sufficiently low-density media such as the one that 
pervades the ISM, kinetic effects such as the Hall effect or 
ambipolar diffusion might also become relevant [ISj . The 
potential relevance of ambipolar diffusion in astrophysi- 
cal dynamos was studied in Refs. [i^l20|. The relevance 
of the Hall effect has been recognized in various astro- 
physical applications [I^, [2ll . [22|, sp ace plasmas |23l - l25| , 
and also laboratory plasmas |26l-|28| . The role of the Hall 
effect on large-scale dynamos subjected to helical forcing 
has also been addressed in the literature [2^, [SO] ■ Less 
attention has received the impact of kinetic effects on 
the small-scale dynamo. A theoretical model of the kine- 
matic small-scale dynamo with Hall effect was presented 
in [3l|, but to the best of our knowledge no numerical 
studies of the non-linear and saturated regime were con- 
sidered in the literature. 

In this paper, we present results from three dimen- 
sional simulations of the Hall-MHD equations subjected 
to random non-helical forcing. The main aim is to study 
the role of the Hall effect in the small-scale dynamo ef- 
ficiency for different values of the Hall parameter. As a 
result of the study, we also discuss the impact of the Hall 
effect on the dynamo saturation values, and on magnetic 
and total dissipation rates. The structure of the paper is 
as follows. A brief introduction to the theoretical frame- 
work known as Hall-MHD is presented in Sect. [Til The 
role of the Hall effect in the efficiency of the dynamo is 



shown in Sect. IIIII In Sect. IIVI we characterize the sta- 
tionary regime that is attained when the dynamo process 
saturates, showing the corresponding energy power spec- 
tra. The energy transfer rates participating in the nonhn- 
ear energy cascade are displayed in Sect. [V] In Sect. IVIl 
we explore the regime of large magnetic Prandtl number 
(i.e., when the viscous dissipation scale is larger than the 
resistive dissipation scale) which, as mentioned, is partic- 
ularly relevant in diffuse media such as the ISM. Finally, 
the conclusions are summarized in Sect. IVIII 



II. HALL-MHD EQUATIONS 

For the sake of simplicity, we consider incompressible 
flows, although compressible effects may be relevant, e.g., 
in the ISM for the formation of structures [32|- Incom- 
pressible Hall-MHD is described by the modified induc- 
tion equation (i.e., with the addition of the Hall current) 
and the equation of motion (the Navier-Stokes equation), 

-TT- = Vx[{U -eVxB)xB]+T]V'^B (1) 

8TJ f B^ 

— = -{U-V)U+{B-V)B-w(p+ — 



F + iyV^U, 



(2) 



where F denotes a solenoidal and non-helical external 
force, which is delta-correlated in time. The velocity U 
and the magnetic field B are expressed in units of a char- 
acteristic speed Uq = \/(U^); f] is the magnetic diffusiv- 
ity, and v is the kinematic viscosity. The parameter e 
measures the relative strength of the Hall effect and can 
be written as 



Ua 



LOpiLo Uo 



(3) 



where Lq is a characteristic length scale, Ua 
•\/(i?^) /iTrnrrii is the Alfven speed 



and 



pi 



yjATTe^n/rrii is the ion plasma frequency (e: electron 
charge, n = Ue = nf. particle, electron, and ion den- 
sity respectively, and m^ : ion mass). Hereafter, we adopt 
C/q = Ua as our typical velocity, thus rendering the Hall 
parameter simply as e = c/(ujpiLo), i.e., a dimensionless 
version of the ion skin depth. 

These equations are complemented by the solenoidal 
conditions for both vector fields, i.e.. 



V-B = 0=: V-C/ 



(4) 



From a theoretical point of view, Hall-MHD cor- 
responds to a two-fluid description of a fully ionized 
plasma: a positively charged ion species of mass rrii mov- 
ing with the velocity field U{r,t), and negatively charged 
massless electrons with the velocity 



Ue 



U-eVxB. 



(5) 



Therefore, from Eqs. ([I])-© we obtain that in the ideal 
limit (i.e. 77 — ?> 0), the magnetic field is frozen to the elec- 
tron flow. As a result, advection, stretching, and folding 
of magnetic field lines (mechanisms relevant for dynamo 
action) are performed by the electron flow rather than 
by the bulk flow, resulting in potential modiflcations to 
magnetic field generation when the Hall effect is not neg- 
ligible. 



III. LINEAR AND NONLINEAR DYNAMO 
EFFICIENCIES 

We performed simulations of the Hall-MHD equations 
with a spatial resolution of 256'^ gridpoints, using a pseu- 
dospectral code |33| . The linear size of our numerical box 
is 27rio (with Lq a unit length), and periodic bound- 
ary conditions in the three cartesian directions are as- 
sumed. We apply the 2/3 dealiasing rule, and there- 
fore the maximum wavenumber resolved by the code is 
kmax = 256/3 w 85. We first consider simulations with 
magnetic Prandtl number equal to unity (i.e., Pm =1). 
The coefficients of viscosity and resistivity in these sim- 
ulations are set to ly = r] = 2 x 10~^, which ensure 
that the dissipation scales are well resolved, i.e., at all 
times the dissipation wavenumbers k^ = (/w^Wi^^)^'^ 
and kjj = (( J^) /rf'Y''^ remain smaller than kmax (here, 
u) = V X [/ is the vorticity, and J = V x B is the current 
density). To evolve the equations in time we use a fully 
explicit second order Runge-Kutta scheme. We note that 
for Hall-MHD, and for velocity and magnetic fields of or- 
der unity, the Courant-Friedrichs-Levy (CFL) condition 
becomes Ai < (Ax)^/e (due to the dispersive nature of 
the whistler waves), which is more restrictive than the 
regular CFL condition Ai < Az. As a result, the time 
step decreases quadratically with the spatial resolution, 
and also decreases linearly with the Hall parameter e. 
The Hall-MHD dynamo simulations are therefore com- 
putationally more expensive than the equivalent MHD 
runs, resulting in the modest spatial resolution consid- 
ered here. 

We first generate stationary hydrodynamic turbulence 
by integrating Eq. ([2]) subjected to random non-helical 
forcing (i.e., such that \7 x F _L F) centered at 
wavenumbers |fc| « fc^ = 3 and delta-correlated in time. 
Once the kinetic energy reaches a stationary regime as a 
result of the balance between the power delivered by the 
external force and viscous dissipation, the hydrodynamic 
simulation is stopped. In a second stage, a random and 
small magnetic field is introduced at small scales, and the 
simulation is restarted with the full Hall-MHD equations 
©-©. 

We performed simulations with different values of the 
Hall parameter e, including a purely MHD case corre- 
sponding to e = 0. Whenever e 7^ 0, a new spatial scale is 
introduced (the Hall scale), which in the spectral domain 
is characterized by fc^ = 1/e. In this paper, we consider 
the cases in which fc^ falls in between the macroscopic 




FIG. 1. Kinetic (thin) and magnetic (thick) energies vs. time 
for e = 0, 0.05, and 0.10 (from top to bottom). 




FIG. 2. Kinetic (thin) and magnetic (tliick) dissipation rates 
vs. time for e = 0, 0.05, and 0.10 (from top to bottom). 



scale kp (set by the external driver) and the dissipation 
scales fcy and A;^, which is the relevant scenario for as- 
trophysical plasmas such as the interstellar medium. In 
such media, the Hall scale is several orders of magnitude 
smaller than the largest scales, and the Hall effect can 
be expected to be relevant only at the smallest dynam- 
ical scales. Note, however, that these arguments ought 
to be regarded as motivations. Although the ordering of 
typical length scales is the correct one, a realistic sepa- 
ration of scales is completely out of reach with present 
computing power. 

In Fig. [T] we show the statistically stationary time se- 
ries for kinetic energy (thin lines) for runs with different 
values of the Hall parameter e. The magnetic energy in 
these runs (thick line) is observed to rise until it satu- 
rates at values which remain a moderate fraction of the 
corresponding kinetic energy. The viscous (thin line) and 
resistive (thick) dissipation rates vs. time are shown in 
Fig. [5] for three runs with different values of the Hall pa- 
rameter. In all these plots, the magnetic dissipation rate 
is observed to grow until it becomes fully comparable 
to the corresponding viscous dissipation rate (even when 
the kinetic energy is larger than the magnetic energy). 

The exponentially fast growth of magnetic energy is 
shown using a lin-log scale in Fig. [3] for the same three 
runs. Note that there is an initial stage where the mag- 
netic field starts growing exponentially fast, regardless of 
the particular value of the Hall parameter e. During this 
early stage of the dynamo, the electron flow is still ap- 
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FIG. 3. Magnetic energy vs. time, showing the exponential 
growth rate during the early linear dynamo regime. 



proximately equal to the ion flow, i.e. Ue ~ U [see Eq. 
dS])]. Keeping in mind that the growing magnetic field 
remains approximately frozen (note that this is strictly 
valid only in the limit 77 — > 0) to the electron velocity 
field, we can anticipate that at some point in time the 
electron and ion flows will start drifting appart from one 
another. Therefore, a second stage arises corresponding 
to a non-linear dynamo (although still "kinematic," in 
the sense that the magnetic field does not affect the bulk 
velocity field), since the magnetic field is being advected 
by the electron flow which at that point becomes a func- 
tion of the magnetic field itself. In Fig. [3] we see that 
although the case e = can be approximated by a linear 
growth rate (indicated by the dotted straight line) all the 
way up to the saturation level, we cannot do the same 
for the cases e = 0.05, and 0.10 since there is a break in 
the corresponding growth rates. This break occurs first 
for the case with larger Hall effect (i.e., e — 0.10), but 
the incremented slope is larger for the case e = 0.05. The 
fact that the dynamo efficiency improves up to a certain 
value of the Hall parameter and then starts decreasing, 
is reminiscent of similar results reported in Ref. |33j for 
large-scale Hall-MHD dynamos. 

To show the relative importance of the Hall term in 
the electron velocity field, in Fig. |4] we display the ratio 
between e\Jk\ and \Uk\ at different labeled times, where 
|Jfe| and \Uk\ are respectively the spectral intensities of 
the current density and of the velocity field at wavenum- 
ber k [note that Ue — U — eJ from Eq. ([5])] . The upper 
frame corresponds to the run with e = 0.05 and the lower 
frame to e = 0.10. The vertical gray line in each frame 
corresponds to the Hall scale k^ = 1/e. In both cases, the 
Hall term becomes gradually non-negligible and eventu- 
ally dominant at the largest wavenumbers of the system. 
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FIG. 4. Ratio between e|Jfe| and \Uk\ at different times (la- 
beled). The top frame corresponds to the run with e = 0.05 
and the bottom frame to e = 0.10. The vertical gray line in 
each frame corresponds to the Hall wavenumber k^ — 1/e. 



i.e., €\Jk\ > \Uk\ at fc > fcg. For the case e — 0.05 (up- 
per panel), the Hall term e|Jfc| becomes comparable to 
|i7fc| at the largest wavenumbers by about t « 24, while 
a similar situation arises for e = 0.10 at t « 18 (lower 
panel). These values of time are remarkably consistent 
with those observed in Fig. [3] for the departure from the 
linear regime in each of the runs. 

In all these cases, there is a third and last stage, which 
corresponds to the saturation of the dynamo. We discuss 
the energy saturation levels in the next section. 



IV. ENERGY SPECTRA 

The distribution of kinetic and magnetic energy among 
spatial scales can be observed in the energy power spec- 
tra Ek vs. k displayed in Fig. [5] for three different values 
of the Hall parameter {Ek is defined such that the total 
energy is E = J dk Ek, and magnetic and kinetic energies 
are such that Ei, + Eu = E) . The thick line in all these 
runs corresponds to the total energy spectrum, which re- 
mains in a roughly statistically stationary regime. Note 
that for these three runs, the kinetic energy spectrum 
remains close to Kolmogorov, i.e., Ek oc k~^/'^, in the 
energy incrtial range. The thin lines display magnetic 
energy spectra at different times, showing the growth of 
magnetic energy as a whole at early times, and satura- 
tion of magnetic field growth at small scales at interme- 
diate and late times. Note that the peak of the mag- 
netic energy in all these cases remains at wavenumbers 
longer than the one where the external force operates 
(i.e., kp — i), as expected for a small-scale dynamo. 

The Kazantsev slope Ek oc k^/"^ Q provides a reason- 
able approximation at small wavenumbers for all these 
cases. Kazantsev's dynamo theory [8| assumes a ran- 
dom velocity field with Gaussian statistics, which is 
spatially homogeneous and isotropic and J-correlated in 
time. Under these assumptions, the two-point magnetic 
field correlation function can be analytically computed 
(see Ref. Q, also Ref. [3J]), and a k^/"^ power law is 
asymptotically expected for the magnetic energy spec- 
trum at the low wavenumber end. Even though Kazant- 
sev's model was devised for pure MHD (no Hall effect) 
and for large magnetic Prandtl numbers, our Hall-MHD 
simulations also reproduce an Ek oc k^''^ magnetic spec- 
trum equally well. This is to be expected, since the Hall 
effect becomes negligible at the lowest wavenumbers (i.e., 
&t k <^ k^ = 1/e). Kazantsev's spectrum has also been 
reported in simulations of small scale MHD dynamos at 
unity Prandtl numbers }ll| . The extension of Kazant- 
sev's model to Hall-MHD in Ref. [31| also recovers this 
spectrum in the regime considered here. 

In summary, a preliminary inspection of the magnetic 
energy power spectra at early times shows no noticeable 
differences between MHD and Hall-MHD. On the one 
hand this is not surprising, since the Hall effect is non- 
linear in the magnetic field, and the magnetic energy at 
early times is much smaller than the kinetic energy at all 
scales. On the other hand, in what follows we show that 
this last conclusion is somewhat premature, since there 
arc other aspects of these turbulent dynamos that clearly 
show the consequences of the Hall effect. 

At saturation, the total magnetic energy reaches a size- 
able fraction of the total kinetic energy, which can be 
estimated within 10% to 20%. More specifically, after 
taking time averages between t = 60 and f = 80 (see 
Figs. [2] and [3]), we obtain the energy ratios E^/E listed 
in Table ID Note that the saturation level of these small- 
scale dynamos, defined as the fraction of magnetic energy 
to total energy in the stationary regime, decreases with 



TABLE I. Global results for runs with different values of the 
Hall parameter e. E is the mean saturation level of the total 
energy, Eb/E is the ratio of magnetic to total energy, kj is the 
average wavenumber for the current density distribution, D is 
the total dissipation rate, and Db/D is the ratio of magnetic 
to total dissipation rate. 
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the Hall parameter. Therefore, although in the linear 
dynamo regime the growth rate increases with the Hall 
parameter e, the magnetic field reaches a smaller satura- 
tion level. 

As mentioned in Sect. IHIl the dynamics of the largest 
wavcnumbers in our simulations is controlled by viscosity 
and electric resistivity. Therefore, the dissipation of mag- 
netic energy mostly takes place in current sheets with a 
thickness which can be expected to be close to the in- 
verse of kri « 85. On the other hand, the width and 
the length of these current sheets will vary from one 
to another [ij, [S^- We can obtain a statistical aver- 
age of the dimensions of our magnetic dissipativc struc- 
tures by computing the power spectrum of the electric 
current density, which is simply k'^Ei,[k). In Fig. [5] we 
show time averaged (between t = 60 and t = 80) cur- 
rent density spectra for three different values of the Hall 
parameter (labeled). All of these spectra are compat- 
ible with a Kazantsev law at low wavenumbers. Note 
that the maximum of these spectra shift toward smaller 
wavenumbers as the Hall parameter increases. Since the 
peak of the spectrum can be associated to an average 
thickness of the current sheets, the above mentioned shift 
can be interpreted as the current sheets becoming rela- 
tively "thicker" as the Hall effect increases. This result is 
in agreement with previous experimental and numerical 
results suggesting that in Hall-MHD the thickness of the 
current sheets is given by the Hall scale rather than by 
the Ohmic dissipativc scale as in the MHD case (see [3g| 
for recent results in support of this interpretation) . 

For each of these runs, we also compute the magnetic 
Taylor wavenumber kj given by 



k", 



Jdkk^ Eb{k) 
Jdk Eb{k) 



(6) 



which arc indicated in Figure [6] by arrows, and arc ob- 
served to remain close, but somewhat to the left of the 
maximum for the corresponding power spectrum. The 
magnetic Taylor scale (i.e., the inverse of kj) can be in- 
terpreted as the mean curvature of the magnetic field 
lines [IJl and of the ensuing current sheets. The value of 
kj also moves towards smaller wavenumbers as the Hall 
scale is increased. The values of kj for each of these runs 
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FIG. 5. Total energy spectrum (thick trace) at t = 72 for 
different values of e (labeled). In each frame magnetic energy 
spectra at t = 18, 36, 72 are also shown (corresponding to 
the thin lines from bottom to top). The Kolmogorov and 
Kazantsev spectra are overlaid (dotted trace) for reference. 



are listed in Table ID In Table U we also list the time 
averaged total dissipation rate D = D„ + Db, clearly 
showing a progressive reduction as the Hall parameter 
is increased. The ratio of magnetic to total dissipation 
Db/D also reduces as e increases, going from approxi- 
mate cquipartition in the MHD case to about 33% for 




FIG. 6. Spectral distribution of current density, i.e., k^Et,{k) 
vs. k for three different values of the Hall parameter (labeled). 
The dotted trace corresponds to the Kazantsev slope fc^". 
The arrows indicate the average wavenumber kj (see Eqn.[B]) 
for each distribution. 



e = 0.10, even though in all these simulations the rela- 
tive content of magnetic energy Eb/ E is comparatively 
much smaller. 



V. ENERGY TRANSFER RATES 

Interpretation of these results on the energy dissipation 
rate requires a detailed analysis of the transfer and con- 

I 



version rate of energy among scales and between the ve- 
locity and magnetic fields, in order to identify the sources 
of small-scale dynamo action in MHD and in Hall-MHD 
turbulence. Equations (II|)-(I21) are knovifn 37| to have 
three ideal invariants: energy, magnetic helicity and hy- 
brid helicity. These are transferred between scales with- 
out losses by the non-linear terms in Eqs. (II])-([2]). In this 
paper we focus our attention in the transfer and conver- 
sion of energy. 



E^\ j d\m^ + \B\')^j 



dk E, 



k 



(7) 



since the non-helical dynamo does not generate heli- 
cal magnetic fields. The dynamo process in this case 
is basically the conversion of mechanical energy into 
magnetic energy by induction, to sustain the magnetic 
fields against Ohmic dissipation. The nonlinear terms in 
Eqs. (HJ-lHl) only redistribute energy (and the other ideal 
invariants) among different spatial scales, in such a way 
to guarantee the global invariance except for dissipative 
losses. To study the energy transfer in Fourier space, we 
follow the method of shell filters originally developed for 
MHD [3^, Ispl ■ which has more recently been extended for 
Hall-MHD [40]. For a given vector field /(r), we define 
/^(r) to be the composition of all wavenumbers between 
K and K + 1, i.e.. 



K-(-l 



/kW 



|fc|=K 



ikr 



(8) 



so that /(r) = E^o/k('')- From Eqs. ^-^ we can 
derive detailed balance equations for the kinetic and mag- 
netic energy in the K-shell 



dtEu{K) 
dtEeiK) 



d^r 



Advection Lorentz 

^ [ -[/. • (£/ • V)Uq + U^ ■ [B- V)Bq 



Dissipation Injection 



fd^r y2\-B^-iU-V}BQ + B^-{B-V)UQ] + rjB ■ V'B, + eJ, ■ {B x Jq) 



(9) 
(10) 



Advection 



Induction 



Dissipation 



Hall 



All cubic terms in Eqs. ([9|)-(fT0|) can be interpreted as 
energy transfer from the third field in the Q-shell to the 
first field in its K-shell, and associated to different physi- 
cal processes as indicated in the equations. For instance. 
Fig. [7^ shows the total shell-to-shell energy transfer rate 
[i.e., the sum of all cubic terms in Eqs. ([5|)- pl]|) ] in the 
steady state of the run with e = 0.05. Light-gray (dark- 
gray) contours correspond to positive (negative) energy 
transfer regions on the (k, Q) plane, located at fractions 



I 

of [0.001,0.010,0.100] of the maximum positive (mini- 
mum negative) value. The first thing to notice is that 
the integral of this function on the (k, Q) plane amounts 
to zero, which is expected to be the case for any con- 
served quantity in the ideal limit. The second aspect to 
notice, is the high degree of concentration around the re- 
gion Q K, Hi, which is indicative of the mostly local nature 
of the direct cascade of total energy in Hall-MHD, just 
as for MHD turbulence (for a detailed study of local and 



non-local contributions to the cascade in MHD and Hall- 
MHD, see Refs. [M El)- In Fig. [71d we show the same 
kind of plot for just the Hall cubic term [sec Eq. (fTO|)]. 
The Hall transfer is non-local, although it is important to 
emphasize that this term is much smaller than the total 
transfer rate. The maximum value for the Hall transfer 
rate is only 4 x 10""*, while the total transfer rate peaks 
at about 0.36. 



The elongated dark spot in the lower part of Fig. [7)3, 
more specifically in the region Q < k^ = 20, indicates 
that energy is transfered backwards by the Hall term 
from small scales to scales larger than the Hall scale. On 
the other hand, the light spot below the diagonal (with 
the dark triangle above the diagonal) for k, Q > fcj = 20, 
indicates that at scales smaller than the Hall scale, the 
Hall term contributes to the direct energy cascade in- 
creasing the transfer rate towards smaller scales. 

We can also compute the energy flux at the wavenum- 
ber |fc| = /c by simply performing 



m) = J2J2 d\\-U^-{U-V)UQ + U^-{B-V)BQ-B^-{U-V)BQ + B^-{B-V)UQ + eJ^-{BxjQ)\ (11 



Each of the five terms on the right-hand side of Eq. (fTTj) 
has a straightforward interpretation. The first term 
(hereafter called H[/[/) is a purely kinetic energy flux, 
which is already present in hydrodynamic turbulence, 
and is responsible of the direct energy cascade in that 
particular case. The second and fourth terms add to zero 
(i.e., WuB + H_B(7 = 0), since they correspond to the ex- 
change between kinetic and magnetic energy by Lorentz 
forces and Faraday induction. The third term (JIbb) is 
flux of magnetic energy associated to the advection of 
magnetic field lines by the velocity field, already present 
in the MHD case. Finally, the fifth term (Hf^") is also a 
flux of magnetic energy, but exclusively due to the non- 
linearity introduced by the Hall current. The first frame 
of Fig. [5] shows the total energy fiux vs. k for the runs 
with e = (black line) and e = 0.05 (gray line). The 
next three frames disaggregate the energy flux into the 
various parts listed above. 

The energy fluxes displayed in Fig. [8] confirm the ob- 
servation made when considering the shell-to-shell trans- 
fer functions. Note that these energy fluxes have been 
normalized by the (time averaged) total dissipation rate 
to allow a comparison between runs with different Hall 
parameter values. The flux associated to the Hall term 
slows down the cascade of magnetic energy toward small 
scales for wavenumbers smaller than the Hall wavenum- 
ber, since H^^" is negative in this range. At the same 
time, the Hall effect speeds up the energy cascade at 
smaller scales, where the H^g' becomes positive. Note 
that the change of sign takes place at the Hall scale (see 
the dashed line in the lowest panel of Fig. ^ . This result 
explains why the current spectrum was observed to peak 
near the Hall scale, and the associated decrease in the 
magnetic energy dissipation rate as the amplitude of the 
Hall term was increased. The slow-down in the trans- 
fer of magnetic energy towards small scales (where it fi- 
nally dissipates) is also responsible for the faster build up 
of magnetic energy through dynamo action in the Hall- 
MHD case. 



r 



Note that although the Hall transfer rate is three or- 
ders of magnitude smaller than the total energy transfer 
rate (see Fig. [7]), the contribution of the Hall term to the 
flux is not negligible. This is the result of the strong non- 
locality of the Hall term: while the Hall transfer is small 
for each value of k and Q, the flux at the wavenumber 
k results from summing over all values of Q, and over 
all values of k < k. The slow decay of the Hall transfer 
far from the diagonal k = Q in Fig. [7] (associated to the 
non-locality) gives a substantial contribution to the total 
flux. 

The other fluxes are also modified by the Hall term. 
The total (i.e., H) and liuu fiuxes decay faster with in- 
creasing wavenumber in the Hall-MHD case, and show 
the build up of (mechanical) flows at scales larger than 
the forcing scale (indicated by the negative value of lijju 
for small wavenumbers). Such an effect for Hall-MHD 
has been predicted from theoretical models [4l| and con- 
firmed by numerical simulations 



Also, the H-uB flux 
peaks at smaller wavenumbers in the Hall-MHD case. 

The modification of the fluxes is consistent with the 
changes in the global dissipation rates. As the Hall term 
reduces the transfer of magnetic energy to smaller scales 
at scales larger than the Hall scale, the global dissipa- 
tion of magnetic energy decreases. Note that this result 
is compatible with studies [42|, |43 , which find faster re- 
connection rates in Hall-MHD simulations with current 
sheets initially set up at small scales. For spatial struc- 
tures such as current sheets, at scales smaller than the 
Hall scale, the Hall term increases the transfer toward 
even smaller scales. As a result, the dissipation rate at 
those particular scales is also increased, which is opposite 
to the result obtained for the global dissipation rate. 



VI. LARGE MAGNETIC PRANDTL NUMBER 

We have so far considered the particular case Pm = 1 , 
i.e., rj = V. However, there are several low-density and 




FIG. 7. Energy transfer rate contour plots on the {k, Q) plane 
(k corresponds to the horizontal axis) for the run with e — 
0.05 at the stationary regime. Both k and Q run from zero 
to kmax = 85. Light-gray filled contours correspond to levels 
at [0.001,0.010,0.100] of the maximum positive value, while 
dark-gray contours display the same levels at negative energy 
transfer rates. The top frame shows the total transfer rate, 
with a peak value of 0.36, and the bottom frame shows the 
Hall transfer rate, with a peak value of 4 x 10~ . 




high-temperature astrophysical plasmas which are char- 
acterized by Pm = v/rj ^ 1, such as the interstellar 
medium, intracluster gas in between of galactic clusters, 
jets, or accretion disks. We performed numerical simu- 
lations with Pm = 10, so that the viscous dissipation 
wavenumber ky becomes much smaller than the resistive 
dissipation wavenumber /c^. As a result, magnetic fluc- 
tuations in this large- Ptti regime may grow at subviscous 
scales. 

In Figure[n]we show the energy spectra for three differ- 
ent runs (labeled). In the top frame we repeat the spectra 
from the purely MHD run with Pm = 1 as a reference. 



FIG. 8. Energy fluxes vs. k, normalized by the total dissipa- 
tion rate D. The black line corresponds to a time average at 
the stationary regime of the run with e = 0, while the gray 
line is for e = 0.05. 



In the central frame we show the spectra for an MHD 
run (i.e., e = 0), but with Pv^ = 10. The thick trace 
corresponds to the total energy spectrum at i = 72, cor- 
responding to the saturation of the dynamo. The various 
thin curves correspond to the magnetic energy spectrum 
at the successive times t = 18, 36, and 72 from bottom 
to top. In the large- Pto regime, the magnetic field grows 
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FIG. 9. Total energy spectrum (thick trace) at t = 72 for 
three different runs (as labeled). Magnetic energy spectra at 
t = 18, 36, 72 (thin lines from bottom to top in each panel) 
are also shown. The Kolmogorov and Kazantsev slops are 
overlaid (dotted trace) for reference. 



in the sub- viscous region of the spectrum. It is apparent 
that for Pm = 10 most of the energy at small scales (i.e., 
large fc's) is magnetic. These results are consistent with 
those reported in Ref. |ll| from very similar non-helical 
simulations. The Kazantsev spectrum Ek « A:'^/^ is also 
a good approximation at large scales, although at large- 




FIG. 10. Magnetic energy vs. time for the three runs as la- 
beled. 



TABLE II. Global results for runs with different values of the 
magnetic Prandl number Pm; t is the amplitude of the Hall 
effect, Eb/E is the ratio of magnetic to total energy, kj is the 
magnetic Taylor wavenumber for the current density distribu- 
tion, and Db/D is the ratio of magnetic to total dissipation. 



e 


Pm 


E 


Eb/E 


kj 


D 


Db/D 


0.00 


1 


0.37 


0.14 


23.1 


0.13 


0.48 


0.00 


10 


0.29 


0.25 


16.4 


0.14 


0.37 


0.05 


10 


0.26 


0.19 


14.3 


0.11 


0.21 



Pm is less surprising, since the kinetic energy is more 
confined toward the small-fc spectral region. 

The lower panel in Fig. [9] corresponds to a simulation 
with moderate Hall value (e = 0.05) and large magnetic 
Prandtl number {Pm = 10). By comparison with the 
case displayed in the central panel (e = and Pm. = 10), 
the dynamo efficiency is larger in the presence of the Hall 
effect, as also confirmed by Fig. [TUl 

Figure [TU] shows the growth of magnetic energy in the 
three simulations in lin-log scale. In the MHD case, the 
linear growth rate at large-Pm (for the same magnetic 
diffusivity rj) is somewhat larger than in the Pm, = 1 case, 
which can be expected as the flow is smoother in the for- 
mer case and turbulent magnetic diffusion is therefore 
less effective. In the Hall-MHD case, the linear regime 
is again found to be followed by a non-linear stage when 
the Hall-effect becomes relevant and the magnetic field 
is advected by the electron velocity, as found in the sim- 
ulations with Pm, = 1. 

Other features of the Hall-MHD small-scale dynamos 
reported in the Pm = 1 simulations can also be identi- 
fied in the Pm, = 10 case. As examples, Table |IT] shows 
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FIG. 11. Spectral distribution of current density, i.e. k^Et{k) 
vs. k for three runs with different values of e and Pm (labeled). 
The dotted trace correspond to the Kazantsev slope k'''^. The 
arrows indicate the average wavenumber kj (see Eqn. O for 
each of the current density distributions. 



saturation values of the total energy, total dissipation, 
and the ratios of magnetic to total energy and magnetic 
to total dissipation for the runs in Figs. |9] and [TOl In the 
MHD case, the increase of the magnetic Prandtl number 
moves the peak of the current density spectrum towards 
smaller wavenumbers (see the values of kj in Table HIl 
and Fig. [TT]). As discussed in Sect. IIV| the HaU effect 
moves this peak further to larger scales. 



VII. CONCLUSIONS 

We present results from three dimensional simulations 
of small-scale dynamo action for magnetic Prandtl num- 



bers Pm = 1 and 10 in conducting flows with the Hall 
effect. This effect is believed to be non- negligible in suf- 
ficiently diffuse media, and its relevance has been rec- 
ognized in various astrophysical, space, and laboratory 
plasmas. As a first step toward a better description of 
dynamo action in such media, only the incompressible 
Hall-MHD equations were solved, and the inclusion of 
compressible effects as well as other kinetic effects such 
as ambipolar difussion is left for future studies. 

However, the inclusion of only the Hall effect acting at 
the smallest relevant dynamical scales of the flow gives 
rise to measurable differences with previous studies of dy- 
namo action. A magnetic non-linear regime is identified 
when the magnetic field (and the current density) be- 
comes large enough to differentiate the electron velocity 
from the bulk flow velocity. After saturation, differences 
in the stationary level of magnetic energy and in the total 
and magnetic energy dissipation rates are obtained, de- 
pending on the amplitude of the Hall effect. Finally, the 
peak of the current density spectrum is found to be de- 
pendent on the strength of the Hall term, with its peak 
moving toward larger scales (smaller wavenumbers) as 
the Hall scale is increased. 

By studying the detailed transfer of energy among 
fields and scales, we observe that the effect of the Hall 
term is twofold: it transfers energy towards larger scales 
for scales larger than the Hall length, and it transfers en- 
ergy towards smaller scales for scales smaller than this 
length. The modification of the energy fiux resulting 
from this transfer is consistent with the observed changes 
in the saturation values of energy and dissipation rate ob- 
served in our simulations. 
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